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Abstract: In this paper we perform a numerical study of the tranverse expansion of 
hadronic scattering amplitudes in the dipole picture of high energy QCD. We go beyond 
the mean field approximation by including fluctuations and also wave function saturation 
effects, and the evolution with both a fixed and a running coupling is investigated. We also 
study the nonperturbative aspects, and as has been predicted earlier, our results indicate 
that the Froissart-Martin bound is saturated once confinement effects are included in the 
evolution. Thus the total cross section increases proportional to the square of the logarithm 
of the cms energy. Using our proton model developed earlier we furthermore see that we 
obtain a reasonable value for the proportionality coefficient. The impact of saturation and 
non-leading effects on this coefficient is also studied. 
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1. Introduction 

The Froissart-Martin (FM) theorem states that hadronic total cross sections must satisfy 



where the coefficient C can be estimated as C ~ 1 /m 2 . The proof of this theorem relies on 
some general properties such as the unitarity of the S-matrix, the existence of a mass gap, 
and the possibility of using subtracted dispersion relations. Although not strictly proven 
within QCD, it is widely believed that this bound should indeed be true for the strong 
interactions. 

The high energy evolution equations (the Balitsky-JIMWLK, or B-JIMWLK hierarchy) 
for the hadronic amplitudes T s {b) were derived in [1], and these equations can also be 
described by the Color Glass Condensate (CGC) formalism [2], or, in the large N c limit, 
by the simpler dipole formalism [3]. The solution to these equations exhibits saturation at 
each b, i.e. T s (b) < 1, with T s (b) = 1 being the black disc limit. This condition is, however, 
not sufficient in order to satisfy the FM bound, since the FM bound is relevant for the 
total cross section which involves an integration over all b. In fact, it is quite obvious that 
if one naively takes the solution to the perturbative evolution equations for T s (i>) , and then 
integrates over all b to get the total cross section as 



one will most certainly violate the FM bound, since the interaction is mediated by massless 
gluons with a Coulomb like behaviour even at large distances. Indeed for an interaction 



crtot(s) < C ■ ln 2 (s/s ) as s — ► 00, 



(1.1) 




(1.2) 
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mediated by a massless particle, such as the photon in QED, the coefficient C in (|1,1| ) is 
infinite, since m 7 = 0. 

In this paper we will study the behaviour of hadronic cross sections with respect to the 
FM bound. We will use our model developed in [4-6], which is based on the QCD dipole 
model, to calculate the growth of the total pp cross section. In [5] we have suggested the 
dipole swing mechanism in order to take into account the missing saturation effects in the 
dipole cascade evolution, and results show that we obtain an almost frame independent 



evolution. As we will later discuss in section |2.2| , the dipole swing has some similarities 
with the saturation mechanism in the CGC formalism. Furthermore, the swing affects the 
expansion of the dipole cascade in the transverse plane, and it is therefore interesting to 
see how large effects it has on the transverse expansion of the scattering amplitude. 

Our main results are presented in section ||. We will see that the confinement mech- 
anism is crucial in order to obtain sensible results. This is especially the case when a 
running coupling is used. Obviously this is to be expected, as otherwise the cascade evo- 
lution favours the formation of too large dipoles. Without confinement the cross section 
grows exponentially in Y (defined as Y = ln(s/so), with so ~1 GeV 2 ), as expected from 
the long ranged nature of the massless gluon fields. 

The leading order cascade evolution is strongly suppressed by the non-leading effects 
[6]. Besides the running coupling, the non-leading effects come from the non-singular terms 
in the gluon splitting function P(z), and the so-called energy scale terms which are related 
to the conservation of p+ and p- respectively [6] (exact energy-momentum conservation 
goes beyond the NLO corrections, however). Once these effects are included (using the 
prescription described in [4]), we see that the growth of the cross section is much reduced. 
Nevertheless, when increasing the energy, one can see that the growth is still faster than 
what is permitted by the FM bound. Interestingly, we will in this case see that atot can be 
fitted rather well by a polynomial in Y, for a s = 0.2 and Y up to around 32 units. 

As the running coupling is a non-leading effect, it might seem strange that the growth 
of the cross section is much faster as compared to the fixed coupling case. This is mainly 
due to two reasons. First, the value of a s gets very large during the evolution (especially 
in the specific model which we use), signaling the breakdown of the perturbative approach. 
Secondly, the total growth is faster due to unrealistically large contributions from larger b. 
However, the growth of T s (b = 0) is actually slower in the running coupling case. 

The transverse expansion of the dipole cascade is of course not really consistent with 
QCD at distances larger than the confinement scale, since it is driven by Coulomb fields. We 
will therefore also study the expansion when confinement effects are modeled by replacing 
the Coulomb propagators, l/k 2 , with screened propagators, l/(fc 2 + M 2 ). In this case a 
ln 2 s growth is obtained, and the FM bound is thus saturated. We also see that we get 



a very sensible result for the coefficient C in (LI). The fact that C can be estimated 
by combining the perturbative growth with nonperturbative initial conditions has been 
proposed in [7] , although the leading order BFKL result gives a way too large value for C. 
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2. The growth of the Black Disc 



Let us consider a hadronic projectile impinging on a hadronic target. The projectile might 
here be an elementary colour dipole, and the target a proton or a nucleus. We will from 
now on denote the generic scattering amplitude between an arbitrary projectile and an 
arbitrary target by Ty(b). When the initial (at Y = 0) projectile or target is fixed we also 
use additional parameters to denote the amplitude. For example if the initial projectile is a 
dipole of size r which scatters off an arbitrary target, then we write Ty(r, b). If we have the 
scattering of two dipoles r and vq we instead write Ty(r,ro,b). The region in b where the 
scattering amplitude Ty{b) satisfies Ty{b) « 1 is called the black disc region, since in this 
case the projectile is strongly absorbed by the target. The region where Ty(b) ~ is on the 
other hand referred to as the "white" region 1 since here the target appears transparent. In 
the region between, the scattering is "grey". Specifically, the black disc region is defined 
as the radius of the disc 2 within which the average amplitude satisfies 

T Y (b) > a for b < R bd (Y), where a « 0.5. (2.1) 

The total cross section can then be estimated as 

a tot (Y) = 2 J d 2 bT Y (b) ~ 2ttR 2 m (Y), (2.2) 

where Rbd(Y) is the radius of the black disc at rapidity Y. (Obviously Rbd{Y) is also 
dependent on the projectile and the target, but we do not explicitely write this dependence.) 
If otot is t° satisfy the FM bound, it is seen that Rbd(Y) can at most grow linearly with 
Y = ln(s/so)- In what follows, we will study the Y dependence for Rm(Y) for different 
situations, using our model developed in [4-6]. 

2.1 The BFKL Growth 

The Y dependence of R^iY) which follows from the solution to the non- linear QCD evo- 
lution equations have been discussed in [7-10], in the context of the Balitsky-Kovchegov 
(BK) equation [1,11] which is a mean field version of the B-JIMWLK hierarchy. A detailed 
numerical study of the BK equation with impact parameter dependence was performed 
in [12]. These works have demonstrated that the BK equation leads to an exponential 
growth of Rm(Y), even if one starts with an initial profile in b which falls off very steeply. 
Using the dipole language, the fast growth of R^diY) can be understood as follows. Assume 
that, for a given dipole projectile of size r, one studies the evolution for b = \b\ >> Rbd(Y), 
where Ty(r, b) << 1. The BK equation can then be replaced by the linear BFKL equation 3 
whose solution in b can be written as 

T Y (r,r ,b) = J d 2 b' J ^n Y (r',b',r )T (r,b\r',b'), (2.3) 

lr The white region can be defined as the region where the scattering between a dipole of arbitrary size 
and the target is small, i.e. where the local saturation scale Q s (b) is smaller than, or the order of, Aqcd [7]. 

2 The average amplitude is isotropic in the transverse plane. On an event-by-event basis, however, there 
is no isotropy. 

3 Strictly speaking one has to be careful when linearizing the BK equation since there might be contri- 
butions to the evolution from inside the black region where the scattering is strong, for details see [7,8]. 
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where ny is the dipole density of the target at Y, and Tq is the basic dipole-dipole scattering 
amplitude. Here we have assumed that the target initially consists of a single dipole of size 
tq. Of course one could imagine more complicated initial conditions for the target, but 
as long as the scattering is weak, the precise choice should not matter. The elementary 
dipole-dipole scattering amplitude, To, can be simplified when the separation of the dipoles 
are large compared to their sizes (remember that we assume b >> Rbd), in which case it 
decays as l/\b — 5'| 4 . It can then be shown that 

T Y (r,r ,b) ~ a 2 s r 2 n Y (r,b,ro) 

- 32 «* Avkn «p k " log— - ) ■ ( 2 - 4 ) 

(irc z Y ) 6 / 2 \ r$r c z Y J 

Here uj = 4ln2-a s , c 2 = 14£(3)a, and a = a s N c /ir. Omitting constants and the slowly 
varying prefactors, we can write 

/ i Qg 2 16fc 2 \ 

r y (r,r ,6)~^^r--|^Lj. (2.5) 

Using definition (|2,1|), it can now easily be seen that this formula implies an exponential 
growth for R^Y), simply because the power- like decay in b, coming from the Coulomb 
fields associated with the exchanged gluons, is too slow to compensate for the fast growth 
in Y. 

2.2 Saturation Effects in The Cascade Evolution 
2.2.1 Multiple Scatterings and Boost Invariance 

The discussion so far has neglected saturation effects in the dipole cascade evolution. In 
the dipole model, unitarization at each b is obtained by taking into account multiple dipole 
interactions. In an eikonal approximation the multiple scatterings can be summed to all 
orders, with the result that T can be written as 

T Y (b) = / 1 - <>xp 7,1 (Xj ■ V> I"' • u ' >) ) • ( 2 - 6 ) 

i j 

Here the brackets denote averaging over different events (notice that it is the event by 
event amplitude which exponentiates). The sums over i and j just denote sums over the 
dipoles in the individual dipole cascades, and the b and Y dependences are implicit in the 
right hand side. Thus we explicitely have Ty(b) < 1 at each b. This condition is absolutely 
necessary for our discussion, as otherwise it would make no sense to talk about the black 
disc limit or the Froissart bound. The Froissart bound determines how rapidly the black 
disc can expand in 6, but if T is not bounded by 1 then there is no black disc limit at all. 

Even though one can unitarize T at each 6, the evolution of the individual dipole 
cascades satisfy the linear BFKL equation. Since what appears as multiple scatterings in 
one frame will appear as saturation effects in the cascade evolution (i.e. in the onium 
wavefunction) in another frame, the formalism is not frame independent. The discussion 
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above should therefore be generalized to the case including saturation effects also in the 
cascade evolutions. Our results presented below includes such effects, and we therefore 
discuss them in this section. 

2.2.2 Dipole Swing 

We have previously argued that one can take into account saturation effects in the dipole 
cascade evolution by including the so-called dipole swing into the formalism. We will here 
once again discuss the idea behind the swing. An analytical proof for boost invariance is, 
however, still lacking. We have therefore so far implemented an approximation in our Monte 
Carlo (MC) code (as have been discussed in [5,6]), and we will at the end of this section 
try to sketch the similarities between our implementation and the saturation mechanism 
in the CGC formalism. 

The generic evolution equations for high energy QCD including all possible pomeron 
interactions are not yet known. Equations have, however, been derived for simple toy 
models which neglect the complicated topology of the full model [13,14]. The saturation 
mechanism in these simpler models have similarities with the CGC formalism, and they 
can also be formulated in a stochastic evolution with similarities to the dipole evolution in 
the presence of the dipole swing, as we now explain. 

These equations can namely be interpreted in terms of (positive definite) k — > k + 1 
vertices [15]. The dipole swing (discussed in more detail below) is a process which instan- 
taneously in Y replaces two initial dipoles, (xi,y±) and (x2,y2), with two final dipoles, 
(xi, 2/2) and (x2,yi). Here (y)x denotes the transverse position of the (anti-)colour end of 
the dipole (here we do not use boldface letters for the vectors). The generic k — > k + 1 
vertices can then be constructed by combining the usual 1 — > 2 dipole splitting with k — 1 
instantaneous dipole swings. In the toy models in [13,14] these vertices give a boost invari- 
ant evolution. However, the amount of information we can extract from the toy models is 
limited. In QCD it is important to take into account the colour degrees of freedom which 
are completely absent in the toy models. The colour structures of the multiple scattering 
diagrams were discussed in detail in [15] and we will here again briefly discuss the colour 
structures. 

The linear cascade evolution is directly related to the leading N c approximation. The 
dipole splitting kernel is proportional to a s iV c = a (for simplicity we here neglect the 
factor -/r in the definition of a as it is completely irrelevant for our discussion), while the 
scattering diagrams are proportional to a 2 . = a 2 /N 2 . To take into account the effects of 
multiple scatterings in all frames one would therefore need to include processes proportional 
to a-a 2n /N 2n (the first factor of a comes from the dipole splitting) in the cascade evolution. 
In the leading N c approximation such processes are absent 4 . To correctly include saturation 
effects in the cascade evolution it is therefore very important that one studies the colour 
structures of the relevant Feynman diagrams. 

4 The reason one includes multiple scatterings is because they are dominant at high energies even if they 
are colour suppressed. The contribution from n pomeron exchange goes like e nujY as compared to single 
pomeron exchange which goes like e" Y , where co is the BFKL intercept. 
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V V 
Figure 1: Two diagrams contributing to interactions within the same cascade as explained in the 
text. The dashed line indicates the cut between the amplitude and the comple conjugate amplitude. 




Figure 2: The respective colour structures of the two diagrams from figure 

Let us consider the situation in figure |l[ Here we have two right moving dipoles, 
(xi,yi) and (x2> 2/2)1 and one left moving dipole (u,v). When rapidity is increased by AY, 
the gluon z can be emitted from (xx,yi) (thus we put the evolution into the right moving 
system). In the left diagram in figure |l| we show just one of the contributions to this 
process. Here z is in the amplitude emitted by the quark located at x%, and is absorbed 
in the complex conjugate amplitude by the antiquark located at y\. In addition to this we 
have two-gluon exchange between the system consisting of xx,yi and z, and (x2i2/2)- In 
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the amplitude a gluon is exchanged between the antiquark at y\ and the quark at X2, while 
in the complex conjugate amplitude a gluon is exchanged between the antiquark at 2/1 and 
the antiquark at 2/2- These processes contribute to the evolution of the wavefunction of 
the right moving system. In addition, there is a two-gluon exchange between the right and 
left moving systems. In the amplitude there is a gluon exchanged between the quark at 
x\ and the oppositely moving antiquark at v, and there is again a gluon exchange between 
xi and v in the complex conjugate amplitude. In the right diagram in figure [j], the gluon 
exchanged between y\ and y2 in the complex conjugate amplitude is instead exchanged 
between x± and yi- 

In figure ^ we show the respective colour structures of the two diagrams from figure |l|. 
Counting the vertices and the colour loops we see that the left diagram is proportional to 
a 5 /iV*. The extra factor l/Nj} comes from the black dots where we project out the colour 
singlet contributions (a dipole is a colour singlet). In the right diagram we can count 2 
loops, and the process shown is therefore proportional to a 5 /N® and thus suppressed as 
compared to the left diagram. 

These diagrams contribute to processess which can be interpreted as follows. The left 
diagram in figure [l] is one of the diagrams which contribute to the process where (xi,2/i) 
first splits into {x\,z) and (yi,z) by the emission of the gluon z, after which the dipoles 
(z, y{) and (x2, 2/2) exchange a gluon, whereby they are replaced by two new dipoles, (z, 2/2) 
and (x2,yi)- Then finally the dipole (x\,z) interacts with the target (u,v). This process 
has therefore a dipolar interpretation, and the step where (z,yi) and (2:2,2/2) are replaced 
by (z,y2) and (x2,yi) precisely describes the dipole swing. The right diagram in figure || 
can on the other hand not be described in terms of dipoles. Here a gluon is exchanged 
between (z,yi) and (2:2,2/2) in the amplitude, but in the complex conjugate amplitude a 
gluon is instead exchanged between (x±,z) and (2:2,2/2)- However, we also see that this 
process is suppressed, since it goes like a 5 /N® instead of a 5 /N%. 

In case the evolution is put into (u,v), the corresponding diagram to the left diagram 
in figure [l] would describe a process where (u, v) splits into (u, z) and (z, v), both of which 
then scatter against (2:1,2/1) and (2:2,2/2)- This is thus a multiple scattering contribution, 
referred to as the "fluctuation" contribution in [16]. If one carefully studies all the possible 
Feynman graphs, then the following picture appears when dipoles in the same cascade are 
allowed to interact. As mentioned above, first one of the dipoles (2:1,2/1) and (2:2,2/2) splits 
(say (2:1,2/1)) into two new dipoles. Then one of the two new dipoles swing with (2:2,2/2), 
and two newer dipoles are formed. Thus we have three new dipoles at the end of the 
process. Any one of these three dipoles can then interact with (u,v). The 2 — > 3 vertex 
which involves the swing goes like a 3 /A^, and is therefore colour suppressed. There are also 
terms which cannot be interpreted in terms of dipole interactions. However, these are all 
proportional to a 3 /N^ and they can therefore be neglected as compared to the dipole swing 
contribution. It was shown in [15] that the colour structures of all the multiple scatterings 
diagrams exactly correspond to the colour structures of diagrams which in the evolution of 
a dipole cascade can be interpreted in terms of the dipole swing and the k — > k + 1 vertices 
mentioned above. Whether or not one can thereby obtain a boost invariant evolution is, 
however, not quite clear. A detailed study trying to attack this problem, and to study the 
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structure of the generated evolution equations is under way [17]. 

Boost invariance (sometimes referred to as the self-duality) of the evolution, and the 
generalization of the B-JIMWLK hierarchy, have previously been discussed in a series of 
papers [16, 18-20]. Note, however, that the k — > k + 1 vertices mentioned above give rise 
to equations more general than the "Pomeron Loop" equations derived in [16, 20] (the toy 
model analogy of this has been discussed for example in [14]). As mentioned, it is, however, 
not know whether or not one can obtain a boost invariant evolution in the full model. For 
example, in [21] it was shown that the higher order corrections arising from the strong 
classical fields in the JIMWLK formalism contains quadrupoles, sextupoles and so on. It 
may therefore be that one needs to include more complicated colour structures in a fully 
consistent formalism. 

In the MC implementation of our model [5,6], each dipole is given one of N% possible 
colour indices, and only dipoles with the same colour index are allowed to swing 5 . Here 
the weight for the swing process was chosen so that it favours the formation of smaller 
dipoles, which explicitely introduces saturation effects in the cascade evolution as smaller 
dipoles both split and interact more weakly. As we discussed in [5] this procedure also 
approximates higher order multipoles, which might be needed in a consistent formulation 
as mentioned above. The similarity with the CGC formalism can then be explained as 
follows. In the CGC, it can be shown that the number density of gluons satisfy [2] 

< - (2.7) 



dY d 2 bd 2 k ~ a s 

due to saturation. If one integrates over rapidity, one rather gets 



iN <X (2.8) 



d'bSk ~ a 2 .' 



S 



What happens in the CGC is that one can continue to pack the hadron with more gluons 
until there are so many gluons overlapping that their mutual interaction is strong enough 
to prevent further occupation at a particular 6. In a semiclassical picture we can think of 
each gluon as a disc of radius ~ l/|fe|- Holding k fixed, one will sooner or later reach a 
point where it is not possible to put in any more gluons of that k. In that case we have 
to increase k, which corresponds to adding smaller discs into the proton. At one stage 
those smaller discs will also fill up the available holes, but there is then more room for 
even smaller discs and so on. One can continue in this way forever, with the typical gluon 
momenta being pushed to higher values, and the total number of gluons therefore never 
ceases to grow 6 . 

The dipole swing works in a very similar way. Since the evolution is driven by the 
1 — > 2 splitting plus the 2^2 swing, the total number of dipoles will continue to grow 
forever. Assume, however, that we wish to put many dipoles of similar size r around the 
same impact parameter b. If the number of dipoles is less than N%, there are no problems 
since the swing is not very likely. However, as soon as we have N% dipoles they can start 



5 The probability that a given colour-anti-colour pair forms a colour singlet is 1/N% . 
6 The production rate of additional gluons does saturate however. 
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to swing, and since in this case they almost sit on top of each other, they will do so as soon 
as the chance is given (the swing favours the formation of smaller dipoles, and in this case 
the swing probability is thus very large [5]). When two dipoles swing they will be replaced 
by two smaller dipoles, with different impact parameters b' . This implies that the dipole 
occupation number satisfies 

dN <JV*~i (2.9) 



d 2 bd 2 r ~ c af 

Here we assume a to be fixed and of order 1/tt, in which case we get a s ~ 1/N C . When 
the number of smaller dipoles around b' gets large enough they will in turn start to swing 
to produce even smaller dipoles and so on. Thus we get a picture which is similar to that 
in the CGC formalism. 

The dipole swing suppresses the growth of the cascade in the transverse plane, and it 
is interesting to see how large effects it has on the evolution. In the next section we will 
therefore present our results with and without the swing. Obviously, we cannot expect 
the swing to change the Y dependence of Rbd(Y) in a qualitative manner, in particular 
we cannot expect it to modify an exponential growth, since the associated interactions are 
still mediated by massless gluons. 

3. Results 

We will in this section present the results obtained from our MC simulation for the growth 
of Rbd(Y). Let us first mention that it is extremely difficult to make analytic predictions 



as in section 2.1. In our simulations we take into account effects of energy-momentum 
conservation, which are related to, but go beyond, the next-to-leading order corrections to 
the cascade evolution. As mentioned in the previous section, we also include saturation 
effects in the evolution, and the full equations which include these effects are not known 
(even in the large N c limit), and once they are known they will probably be extremely 
complicated to solve. Furthermore, analytic estimates can never fully take into account 
the full impact parameter dependence of the evolution, which is relevant for the present 
study. 

The first results we show are for the case of a running coupling, and without any 
confinement effects. In order to avoid singularities, the value of a s (p±) is frozen below 
the scale p± = 2/r max where r max is a free parameter which in our full model sets the 
confinement scale. In our previous studies [6] we have set r max = 3.5 GeV -1 which is also 
the value we will use in the present study. (Aqcd is fixed to 0.22 GeV) 



In the discussion in section 2.1, we considered the projectile to be an elementary 
dipole. Due to the fact that higher Y values are extremely time consuming to simulate, 
we will here use a projectile which is more dense initially. We will therefore consider 
pp scattering where the initial proton is modeled as consisting of three dipoles (with a 
Gaussian distribution in sizes determined by the scale r max ) in a triangular configuration 
as was discussed in [5,6], where the frame independence of the process have also been 
demonstrated. We will therefore compute Tyib) in the CM frame, since this is numerically 
the least time consuming frame. 
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Figure 3: The transverse profile of the scattering amplitude T including a running coupling, and 
without confinement effects. The left figure excludes the swing while the right figure includes it. 
In both figures the lowest curve is calculated at Y = 8, and Y is increased by 2 units for each new 
curve. In both cases an exponential growth of Rbd(Y) is seen. The main mechanism driving the 
growth is clearly the very fast growth of the white region. 
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Figure 4: The black disc radius Rbd(Y) plotted as a function of Y, as calculated from figure 
|^. The squares correspond to evolution including (lower set), and excluding (upper set) the swing. 
Together with each curve, the corresponding fits are also shown. 



Our first result are shown in figure ||. Here we plot Ty(b) as a function of b (we 
average over the angle). We show two plots, including (right plot) and excluding (left 
plot) the swing. The lowest curves are calculated at Y = 8, and Y is increased by 2 units 
for each new curve. We see very large contributions from large impact parameters, and 
the resulting profile is very flat. In this case Rbd(Y) grows exponentially which can be 
seen in figure |4j. Here, the upper curve excludes the swing, and we find that Rbd(Y) can 
be fitted as R b d(Y) « 1.8-exp(0.15y)GeV _1 . If the dipole swing is included we instead 
get Rbd(Y) ~ 2.0-exp(0.13y)GeV~ 1 . Also in the case where we include the swing, we see 
that RbdiY) grows very rapidly. As mentioned in the introduction this is partly due to the 
breakdown of perturbation theory since the value a s (p±) is very large during the evolution, 
which is moreover extremely sensitive to the infrared cutoff r max . 

Of course, since in this case the initial dipoles are rather large (with sizes close to 
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Figure 5: The transverse profile of T for onium-onium scattering, without the swing. We again 
find a very flat profile, with very large contributions from large b. The lowest curve is calculated at 
Y = 12, and Y is increased by 4 units for each new curve. 

r max) one is in the soft region already from the beginning. We therefore also study onium- 
onium scattering, where the initial onia have small sizes, we take the two cases ro = 0.5 
and 1 GeV -1 respectively (both initial onia have the same size). In figure [| we show the 
transverse profile for the case ro = lGeV -1 . Here we do not include the swing and we once 
again find a very flat distribution implying an exponential growth for R^diY). The result 
for the case ro = 0.5GeV _1 looks essentially the same. 

Next we switch to a fixed coupling, a s = 0.2. It is well known that the leading 
order dipole cascade also in this case shows a fast diffusion towards large dipole sizes. 
In [4] we demonstrated the very large effects of energy-momentum conservation on the 
evolution. In this case the production of both small (from p + conservation) and large 
(from p_ conservation) dipoles are suppressed, and the growth of the cross section is 
severely dampened. 

The transverse profile of the scattering amplitude is shown in figure ||. As compared to 
figure H we see that the growth is significantly reduced. This is because the unrealistically 
large contributions from larger b present in the running coupling case are much reduced. 
From figure ^ it seems that, for rapidities under Y ~ 20, the growth is even linear 7 . 
However, as Y increases, we clearly see a deviation from the linear growth, which becomes 
greater for the largest Y. Clearly, the FM bound-breaking growth arises not because the 
central regions reach the black disc limit too fast, but rather because the white region, 
where the scattering is very weak initially, turns grey too fast. This indeed confirms the 
expectation that the FM bound is in this case violated due to long range contributions 
coming from the perturbative Coulomb fields [8]. What is interesting, however, is the fact 
that the growth of Rbd(Y) is much slower than expected from the leading order evolution. 
This is clearly illustrated in figure |. What we see here is that the shape of Rbd(Y) cannot be 

7 In order to see the FM bound-breaking growth sooner, we have here defined Rbd for a = 0.3 in ( |2.l[ ). 
This will of course affect the absolute value of Rbd, but it does not change the fact that Rbd grows more 
than linearly with Y. 



-11- 





2 4 6 8 10 12 14 
b(GeV" 1 ) 

Figure 6: The transverse profile of the scattering amplitude for a fixed coupling, a s = 0.2. In 
this case the swing is not included. The lowest curve is calculated at Y = 8, and Y is increased by 
2 units for each successive curve. Although the growth of Rbd(Y) seems to be linear for lower Y, 
we can see that the white region grows very fast, implying an increasingly faster growth of Rm{Y) 
as Y increases. 



> 

CD 

O 

.O 

cc 




Figure 7: The black disc radius Rbd(Y) as a function of Y as calculated from figure ^ (squares). 
In this case Rm(Y) can be fitted by a polynomial (the dashed line) as explained in the text. 



fitted by an exponential, at least for Y up to 32. The growth accelerates as Y increases, so as 
Y — > oo the growth should eventually reach an exponential (indeed the white region already 
grows exponentially). For the present values we instead find that the shape can be fitted by 
a polynomial, and in figure @ we show a fit R = (-3.3+0.76-y-0.02-y 2 +6-I0~ 4 -y 3 )GeV^ 1 , 
which is the best fit we have found. Such a fit would imply a cross section growing like 
a tot ~ln 6 (s/so)- This growth is further reduced if we include the dipole swing, as we 
illustrate in figures ^ and ^. Even if the growth seems to be linear, one can again see that 
the white region expands rapidly and the growth of Rm accelerates as Y increases. The fit 
in figure |9| corresponds to a quadratic fit. 

In order to satisfy the FM bound we must introduce a scale above which the gluon fields 
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Figure 8: The transverse profile for fixed coupling and including the swing. The lowest curve is 
calculated at Y — 10 and Y is increased 2 units for each new curve. 
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Figure 9: The black disc radius RbdiX) as a function of Y as calculated from figure || 

falls off exponentially. In [6] we replaced the momentum space Coulomb propagators, 1/A: 2 , 
by screened propagators, l/(fc 2 + M 2 ), where M = l/r max . The dipole splitting kernel, for 
the process (x,y) —> (x,z) + (z,y), is then modified as 



{x -yf 



(x - z) 2 (z - yy 



1 x - z 



T*m.n.T. X Z 



Ki(\x-z\/ 

i max ) 



1 z-y 

max \% y\ 



K 1 (\z-y\/r r> 



(3.1) 



Similarly the dipole-dipole scattering amplitude To(x.y\u. v) is modified as 



ai, o [ \x — v\\y — it 
— -In 



\x — u\\y — v\ 



2 



K Q (\x -u\/r max ) - K (\x -v\/r % 



max )+K (\y-v\/r max ) 



(3.2) 
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Figure 10: The transverse profile of T including confinement effects via equations (3.1) and ( |3.2| ). 
In the left figure the swing is excluded while it is included in the right figure. The lowest curves are 
calculated at Y = 8, and Y is increased by 2 units for each new curve. Here the growth is clearly 
linear, and as expected we see that the growth of the white region is considerably suppressed. 



Here Kq and K\ are modified Bessel functions which both behave as K (x) ~ y \ e ~ x f° r 
large x. In the analysis of section [2.1| , this would imply an exponentially decaying profile 
in £>, which would compensate the exponential growth in Y of the BFKL solution. The 
evolution should then satisfy the FM bound. This has been used in [7] to estimate the 
constant C in Ql.lQ as 

C = 2"(£f, < 3 ' 3 > 

where u> was defined in ( |2.4j ), and fi is the confinement scale, i.e. r^ ax in our model. The 
constant C can thus roughly be determined by combining the hard pomeron intercept with 
the non-perturbative confinement scale, albeit in a heuristic fashion. 

Recently it has been shown that experimental results on pp collisions favours a ln 2 s fit 
(rather than a Ins fit) to the total pp cross section [22]. The fit in [22] has the form 

^co + Cl ln(^) +C2 ln 2 (^), (3.4) 

where we have neglected terms which fall off as a power of s. The various coefficients 
above were found to be cq w 37 mb, c\ pa —1.4 mb and c 2 ~ 0.28 mb. In [7], the value of \i 
has been argued to be around 2m n ~ 0.28 GeV, which interestingly is equal to with 
r max = 3.5 GeV -1 . In this case the value of ( |3.3; ) would be around 9 mb, considerably 
higher than c 2 above. Of course, the value of ( |3.3|) is not supposed to reproduce c 2 since it 
has been derived under rather crude assumptions (any attempt to include nonperturbative 
effects will admittedly be heuristic as well). It is also a well known fact that the leading 
order BFKL exponent uj is too large to fit data. If we would for example replace oj by its 
NLO value, oj ~ 0.3, C would be reduced almost by a factor of 4. 

In figure |To| we show the transverse profile of the scattering amplitude, calculated 
using ( jO| ) and (|3.2[) , and using a running coupling. Note that, as compared to figure |6|, 
the difference now is that the growth of the white region is considerably suppressed, as 
expected from confinement. In this case we can see a linear growth, which is illustrated in 
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Figure 11: The black disc radius Rbd(Y) calculated from figure 10 excluding (upper set of squares) 
and including (lower set of squares) the swing. Lines represent linear fits to the squares. 



figure 11. Here we find that Rf,d can be fitted as R bd (Y) = (-1.36 + 0.39y)GeV" 1 when 
the swing is excluded, and as Rbd(Y) = (—1.70 + 0.3iy)GeV~ 1 when it is included. This 
would imply a cross section growing as a ~ 2it • 0.39 2 ln 2 s GeV 2 = 0.37dn 2 s mb for the 
former case, while for the latter case we would have a ~ 2vr • 0.31 2 ln 2 s GeV 2 = 0.24dn 2 s 
mb. The results are indeed quite close to the C2 term in ( |3.4[) . We might also ask how good 



the approximation a ~ 2itR 2 d (Y) is. The MC results for the transverse profiles in figure 



10 can be estimated rather well (for fixed Y ) by simple Gaussians, T = T(0) exp(— c • b 



Thus, approximating these curves by such functions the cross section would be given by 

ai = 2 J d 2 bT(0) exp(-c • b 2 ) = ^T(0). (3.5) 



On the other hand, Rbd is defined as T(Rbd) = a with a = O(0.5), and the approximation 
a ~ 27r Rl d {Y) then gives 

^=^jm\ (3. 6) 

c \ a J 

Thus 

- * jm\ (3 . 7 ) 

cji T(0) \ a J 



For the Tevatron for example, we have T(0) ~ 0.7, while in figure |ll| we have set a ~ 0.6. 



This would give 02jo\ ~ 0.5. For higher energies where T(0) -> 1 we get a"2/ci ~ 1 so 
that the approximation a ~ 2irR 2 d (Y) works reasonably well. 

Of course we might as well directly calculate the total cross section using the MC, such 
as we did in [6] . Thus we can try to fit a curve of the form (|3.4p to our results [6] for the 



8 Given the fact the functions K\ in the kernel ( |3.l[ ) fall exponentially one might think that an exponential 
function would better approximate the 6-profile. The good Gaussian approximation is related to our initial 
proton model which has a Gaussian distribution [6] . Actually our fo-profile has a somewhat longer tail than 
a Gaussian for large b due to the fluctuations in the evolution, as was discussed in [6]. 
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total pp cross section. We thus parameterize the cross section as 

a = A + B\ns + C\n 2 s, (3.8) 

where s is measured in units of GeV 2 . In this case we find the results A = 34 mb, B = —1.8 
mb and C = 0.30 mb which are quite close to the values cq, c\ and C2 in (|Q|). We also 
try a fit a = D + E ln 2 s which works equally well, and we find D = 20 mb and E = 0.24 
mb. We thus find that the Froissart bound is saturated, and by combining the perturbative 
evolution with nonperturbative confinement effects we moreover obtain a value for C which 
is consistent with data. 

4. Conclusions 

We have in this paper studied the growth of the black disc radius Rm{Y) for hadronic 
collisions, using our model developed in [4-6]. Using a purely perturbative approach one 
cannot expect the FM bound to be satisfied, and for the running coupling case we indeed 
find an exponential growth for R^Y). On the other hand when a fixed coupling is used we 
again find a fast growth of R^iY), but the growth is in this case not an exponential as in 
the running coupling case, but can be fitted by a polynomial (at least for rapidities up to 32 
units, corresponding to s ~ 10 14 GeV 2 ). This is due to our inclusion of energy-momentum 
conservation effects in the evolution, which severely dampens the leading order growth. 
However, it can be seen that the white region expands exponentially, which should imply 
an exponential growth for Rm(Y) eventually. The fact that the running coupling case 
shows such a fast growth is because, without any suppression of large dipoles, the coupling 
gets very large during the evolution. Moreover, one gets unrealistically large contributions 
from large transverse separations. 

We model confinement effects by replacing the Coulomb propagators by screened prop- 
agators, in which case R^diY) grows linearly with Y, implying that the FM bound is actu- 
ally saturated. Furthermore, including saturation effects during the cascade evolution, we 
see that we obtain a value for the coefficient C in ( |1.1| ) which is consistent with data. 

Obviously, our specific model used here makes sense only when including nonpertur- 
bative effects since the initial dipoles in the proton are quite large. However, also in the 
case where we start from a, smaller, single dipole did we see a very fast growth. It is also 
interesting to see what would happen if one starts with a system containing perturbative 
dipoles, which is at the same time quite dense. Then we would expect saturation to slow 
down the evolution, although we will of course still get an exponential growth if no confine- 
ment effects are included. An interesting initial model is a proton consisting of three "hot 
spots", i.e. three saturated spots inside the proton. This model was discussed in [9]. The 
sizes of these spots, related to the scale of chiral symmetry breaking, is estimated to be 
around 0.3 fm. Thus one may at least initially neglect confinement effects for the evolution 
of each spot. 

To test the sensitivity of our model to the initial assumptions, we have also tried a 
model where the proton initially consists of 6 dipoles in 3 spots (2 dipoles in each spot), 
where each spot has a size around 0.3 fm. Including the swing, and with a running coupling 
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and no confinement effects, we again find a very fast growth of atot- For spot sizes of around 
0.3 fm, we find that the cross section shows a s ' 21 behaviour. 

Finally, when confinement is included we find that a fit D + £Mn 2 s gives, E ~ 0.31mb 
which is higher than the result 0.24 mb found above. However, we should mention that 
this initial model also reproduces the available high energy data rather well. Furthermore 
we also get good results, as in [6], for the diffractive and elastic cross sections, both in 
pp collisions and in DIS. To the accuracy of our model, we therefore get a reasonable 
description of data also with such an initial model. 
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